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We study the quantum correlations of the radiation emitted by three level atoms (cascade type) 
interacting with two driving fields. In the linear regime, and in the Weisskopf-Wigner approximation, 
we show that the atomic and the two-photon density matrix are equivalent to each other. This 
facilitates the tomography of the two mode state to be realized by measurements on either the atomic 
system or the emitted fields. While, in general, one needs 4^^ measurements for the tomography of 
a A'^ photon state, we show that one needs (A^ + 1)^ — 1 observables for the tomography of photons 
' emitted by an atomic system. Thus there is an exponential reduction in the number of observables 

T-H ' for the reconstruction of the class of A'' photon states emitted by atoms. We show that the driving 

. field strengths and detunings provide the control parameters for the preparation of a specific target 

04 ' state. Finally, we study the characterization of entanglement of the two photon state. We observe 

^ ' that a characterization of entanglement in terms of a single parameter is not possible when the 

^ , system is in a mixed state; therefore, we provide a description in terms of the newly introduced 

■ probability distribution for entanglement, in various regimes of interest. 

G> ' PACS numbers: 

(N : 

^! I. INTRODUCTION 

> ■ 

d ' Quantum state tomography (QST) [l[ is the first step in any quantum information processing (QIP). For a photon 
^ \ qubit, the state determination may be in any of the degrees of freedom, eg. polarization, phase space distribution, 
O^. or occupancy in the Fock space At the experimental end, reconstruction of the quantum state is either by 
direct photon counting measurements or indirectly by optical homodyne /heterodyne measurements - as in the 

ps) . case of continuous variable tomography. James et al Q have proposed a general scheme for the reconstruction of 
^ ■ the polarization state of a pair of entangled photons generated in type-II parametric down-conversion (PDC). The 

On , PDC two-photon probability distribution has also been determined experimentally 0] by measuring the joint photon 

■ probability distribution. Theoretically, Vogel and Risken Q have shown that the Wigner function can be determined 
by optical homodyne detection (OHD), which was subsequently followed by the experimental realization of the same 
0. The phase-space quasi probability function has also been determined @ by photon counting measurements. Single 

] photon [3 and two photon Fock state tomography have been experimentally studied by Lvovsky et al Q • Spin state 
tomography has been studied using group theoretic methods [9|. So far, tomography studies on photons have been 
mainly in PDC systems. 

Preparation of a specific target state as required by the protocol of interest is also necessary for any QIP. Control 
of atomic sources for the realization of specific photon states may be achieved via coherent atom photon interactions 
where the driving field strengths and the detunings of the applied fields are the control parameters . An inevitable 
; ^ > consequence of atom photon interactions (coherent) is the existence of strong quantum correlations [lol - [l2j . Coherent 
I driving fields interacting with the atoms give rise to strong atomic coherence which then gets transferred to the 
emitted radiation. Quantum correlations, particularly the entanglement aspects, have been the most studied but not 
easily understood, especially when the state is not pure. On the other hand, the preparation of pure entangled states 
itself poses a challenge. While PDC sources offer good fidelity the conversion efficiency is, however, very low since it is 
a second order process. Further, the input field strength is the only control parameter. A further disadvantage is that 
there is no control over the emission statistics. On the other hand it is well known that the statistics of emission in 
resonance fluorescence can be made sub-poissonian or super-poissonian. In other words, there can be programmable 
delays between successive emissions. The disadvantage of monitoring emission in resonance fluorescence is that the 
emitted radiation is over 47r and the signal to noise ratio per solid angle may be small. This may be overcome by 
coupling to a cavity which is tuned to the mode frequency of emission. 



o 
o 



X 



*EIectronic address; lsns@iitk. acini 
t Electronic address: Vravi@iitk.ac.iii1 



2 



In this paper we consider three level cascade systems interacting with two driving fields. We study the quantum 
correlations of the radiation emitted by these systems in free space. In particular, we address three issues : 
(i) quantum state tomography of the two-mode emitted radiation: In the far-field approximation (lowest order in l/r, 
and t > r/c) and using the Weisskopf-Wigner theory of spontaneous emission, we show that it is sufficient to measure 
a set of 8 observables to completely determine the two-photon density matrix . This is because a iV - photon state of 
the emitted radiation in our system is restricted to a Hilbert space of dimension {N + 1)^ — 1 instead of the full Hilbert 
space of dimension 2-^^. Accordingly, there is an exponential decrease in the number of observables that are needed for 
tomography. For the present case, iV = 2, and the number of required observables is 8. An additional feature is that 
there is a one-one mapping of the atomic and the photonic density matrix and therefore a knowledge of one leads to 
the other. Unlike conventional methods of tomography where measurements are restricted to the emitted radiation 
only, in our scheme there is no need to perform the measurements entirely on either the emitted radiation or on the 
atomic system. We identify six observables of the atomic system and two observables of the emitted radiation as a 
convenient set. Thus this approach has an edge over the conventional methods of tomography. 

(a) Control and preparation of specific states: As already mentioned, it is not easy to prepare entangled systems with 
a high degree of purity, and with high efficiency. We show that three level cascade systems provide a very good source 
of entangled two mode emissions with very large probabilities. In fact, even in the steady state one has the option 
of preparing almost pure fiducial states. Bell states with reasonable fidelity, or mixed states by suitably setting the 
driving field strengths and detunings. 

(Hi) Characterization of mixed state entanglement: Conventional entanglement measures such as concurrence and 
negativity are known to provide an incomplete description of entanglement in the mixed state regime. We evaluate 
the probability density of entanglement for the two-mode state following Bhardwaj and Ravishankar [l3| . 

The plan of the paper is as follows: In section II we describe the model, and set up the equation of motion. In 
section III we dcmonstarte the reconstructibility of the atomic density matrix and also the photon density matrix 
by a suitable set of measurements. In section IV we discuss the control of pure/mixed state preparation and the 
characterization of entanglement. 



II. MODEL 



Emission from three level cascade sytems are of interest, especially for the strong quantum correlations that they 
exhibit. The earliest classic work of Clauser llSl in Cs atoms demonstrated the non-classicality of the emitted radiation 
in these sytems. The resonance fiuorescence [l6[ and the emission statistics 17[ in this system have been well studied 
and the behavior is qualitatively predictable in most regions of the parameter space. In view of t his, and the large 
number of recent experiments on electromagnetically induced transparency (EIT) in vapours [l^l, we choose to 
study three level cascade systems. In passing, we would like to mention that the analysis presented in this paper is 
equally applicable to 'A'and 'V systems by properly taking care of the detunings. The scheme, shown in Fig.l, has 
only two allowed dipole transitions of energy separations uji and uj2 corresponding to the |2) — >• |1) and |3) — > |2) 
transitions. Two counter propagating ( Doppler free geometry) driving fields of nearly equal frequencies ujli and u!l2 
and respective strengths Jli and ^2 are resonant with these two transitions. The decay constants of the energy levels 
1 3) and |2} arc indicated by and T2 respectively. The parameters Ai, A2 refer to the detunings of the driving 
fields. 



A. Equations of motion 



Since the input parameters are known - and so is the dynamics - it is straightforward to solve for the state of the 
atomic system. There is, however, no corresponding procedure (availability of a master equation) that determines the 
state of the emitted radiation. One way is to determine the atom-|-radiation density matrix, and then trace over the 
atomic part which is cumbersome. In this paper we will show that this tedious procedure can be circumvented. We 
further show that the atomic density matrix determines the photon density matrix completely! 

We start by making an ansatz for the atom+field state as 



I* W) - E Wl*! 7); * = 1, 3; 7 = 0, 1, 2, 3; (1) 

where i denotes the atomic index while 7 denotes the photon mode index. We alternately use the notation |i; 7) = 
\i; 72, 71). The notation \i] 7); 7 = 0, 1, 2, 3 is more commonly used in quantum information where 7 ~ 7271 in base 2. 
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FIG. 1: (Color online) Three level cascade system corresponding to the ^"^ Rb atoms driven by two fields lji and UJ2- 



More explicitly, 7 = 2°7i + 2^72. In writing the above equation we make the following assumptions whieh are eentral 
to our analysis: (i) The rotating wave approximation (RWA) which restricts the transitions sueh that a photon is 
created only when an atomic excitation is annihilated and vice- versa. This is very similar to the ansatz of Weisskopf 
-Wigner theory of spontaneous emission, (ii) The single-photon approximation |14l | - the rise time or the excitation 
time of the atom is inversely proportional to the driving field strength. When D, is small compared to the atomic 
decay constant F, the rise time is large and the atom completes fewer cycles of excitation/de-excitation in an 
atomic lifetime. This situation enables us to assume that for small times {t < atomic lifetime) the Fock space of the 
two modes of radiation has at most one photon in each mode. Thus 71, 72 = 0, 1 only. We therefore restrict our study 
to the situation where the ground state excitation is weak. 

The interaction Hamiltonian in the interaction picture and in the RWA and the dipole approximation has the form 



h.c.) + hVt2[e~ 



h.c.) 



(2) 



Here we treat the driving fields classically and Hi 



E, 



1,2, are the Rabi frequencies corresponding to 



the two driving field strengths, = \i + 1) = \i) (i + 1|; i = 1,2 and af = |i + 1) (i + 1| — \i) {i\ arc the 

atomic transition operators. The equation of motion for the 12 complex amplitudes aj is then given by 

^na]{t) = {^■^\H\^{t)) 

It might appear that, in general, the two-photon density matrix elements corresponding to the two modes 



(3) 



F 
"11' 



(4) 



would require a set of 16 observables for its full determination. However, by employing an equivalent approach in the 
next section, we show that not all of the complex amplitudes al{t) are independent since they are constrained by the 
production process and selection rules. We further show that a\,a\, a^, a|, aj, are zero because of RWA . 



4 



III. RECONSTRUCTION OF THE TWO-PHOTON STATE 
A. Equivalence of the atomic and photon density matrix 

In general, the field E{r, t) emitted by the atom has the mode expansion 

3+Hr,t) = i^ffcefc^Aflfc^AWe'^*'-'- (5) 

fc,A 

for the positive frequency part where £k = {2'!Thck/Vy^'^ for a given mode k and A is the polarization index. We 
now assume that the fields emitted by the two transitions to be dominantly centered around fci and fc2 respectively 
and are defined as in eq(2.12) of Titulaer and Glauber [l^. The creation and annihilation operators satisfy the usual 
commutation relations [ak,\{t),al., yit)] = Sk,k'^\,\' for a given mode k. In order to establish the equivalence of the 
atomic state with that of the radiation we consider the auxiliary state in the full Hilbert space of atom+radiation 
defined by 

= EE E B'pip',q'i4^t)r{al{t)n;0k,;^ 

i i' p,q,p' .q' 

(6) 

In the single photon approximation, each mode has either no photons or at most one photon, i.e, {p, q,p', q'} = {0, 1}. 
The indices ki,k2 correspond to the two modes and \0k2',0ki) is the corresponding vacuum. Recall that p^^^{t — 
Q) ^ ® . Also recall that the annihilation operators for the two different modes are mutually commuting at a 
given time, i.e., [a^^ (t), a^j (<)] = 5kik2 8'iid likewise for the creation operators. The two-photon density matrix may 
now be obtained by taking a partial trace over the atomic variables to be 

p^(t) ^ ^ < >= ^P,?;P'..'(<W)"(«l(i))'|0^.;0^2> {Qk,M{ak2it)Y' {akAtW' 

(7) 

where the coefficients Bp q.pi qt = q p' q'- This expression is equivalent to the form that would follow from eq.l, 

as may be seen by observing that |i;7) = |i;7i,72) — (aj)'''^ (aj)'''^ |i; O2; Oi). Note that this is very similar to the 
definition of the most general density operator of n photons in a single mode given by eq.2.15 in [l9[. Following the 
Weisskopf-Wigncr theory of spontaneous emission, we stipulate the far field approximation wherein the field operators 
i?+(r, t) of the emitted radiation at the point of the detection r arc proportional to the atomic spin operators at the 
retarded time aj (t — |r — ro|/c). In the far field for an atom located at Tq, terms proportional toO(l/|r — rop) may 
be neglected [22, l23l- Thus, for a given mode k of the source field, 

i?r)(r,t) . E-ir,t) + ^^^J^^^kit - |r - r,\/c) (8) 

where i?(^(r,i) is the vacuum field contribution and ry is the detector efficiency. In this paper we consider the ideal 
situation of 77 = 1. The first term due to the vacuum does not contribute to any signal detection and may be ignored 
[2l| for detection times larger than the inverse of the optical frequencies. We use the notation = dkk+i and 
= Ok+ik- For convenience we replace the mode labels fci and k2 with '1' and '2' respectively. It is convenient to 
introduce the auxiliary state 

fit) = 

^ /2(r2)/i(ri)^ ^ i?p,,,p-,,,(a2+(i2/c)r(a+(ti/c))«|z;02;0i) (z;02;0i|(ar(ti/c))''(a2-(i2/c)r' (9) 
i p^q.p' ,q' 

2 I 

where fk{rk) = i~ ^TieoSk ^ ~ ^' ^^"^ ~ i^Ti/c]i = li 2. The atomic operators satisfy the commutation 

relations [aij,<Jki\ = Sjkcru — Suakj which follow from the identity (jij(t)uki{t) = crii(t)Sjk for the atomic operators 
CTy (0) = \i) When the point of observation is the same for both the modes r2 = ri = r. Therefore 

f{t)^f2{r)h{r)Y^ ^ Bp,qy^q,{a+{t-r/c)r{a+{t-r/c)n;02;0i) (z; O2; OiKa^ (i-r/c))«' (a-(t-r/c)f ' (10) 
i p,q.p',q' 
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Again since the atom and photon density matrices are separable at t=0, the summation term in the r.h.s is nothing 
but the atomic density matrix. Thus, it follows that 

p'{t)®T^ ^Npl®p^{t~r/c), (11) 

where M = f2{T)fi{r) is the geometric factor and = |0i; O2) (Oi; O2I. The normalized density matrices are obtained 
by dropping N . 

Eq. 11 merits some more explanation. First of all, p'^[t) at any time gets determined completely by p^{t — ^), the 
atomic state at an earlier time t — r / c. More interestingly, we observe that p'^ is of rank three. This reduction in the 
rank is a manifestation of Schmidt decomposition as applied to the atom+ radiation system. More specifically, the 
dynamics - which leads to the equivalence of the atomic (retarded time) and the field operators, yields the following 
important relations: (i) p'^ belongs to the space orthogonal to the state |2)^ = |10). (ii) Let us set up matrix elements 
for p^ in the basis {|00), |01), |11)}. Eq. 11 asserts that the matrix in this basis is identically equal to the matrix 
for the atomic state set up in its standard basis labeled by the energy eigenstates {|1), |2), |3)}. Explicitly, the 
normalized photon density matrix has the form 



1 / E.a°«r E,«°«f \ 

p'{t) = T7 E,«.'«r nM? E,«.'«r 

( Pll 9^2 Pts \ 

^ Pi2 P^3 =p\t~rlc) (12) 

where the subscript emphasizes that the atomic density matrix is evaluated at a retarded time t,. = t — r / c. The 
atomic density matrix may be determined by solving the Liouville equation 



ihp^[H,p]-'-^{V,p} (13) 
where H is the Hamiltonian given by eq(2) and T is the relaxation matrix. The second term is explicitly given by 



B. Measurements and Observables 



One of the important consequences of the equivalence of the atomic and the photon density matrix, shown in Eq. 
12, is that it gives us the freedom of choosing measurements on either the atomic or the photons emitted depending 
on the practical feasibility. The existence of more than one complete set of measurements further helps in minimizing 
uncertainties in the determination of the quantum state. The complex amplitudes, and hence the density matrix 
elements, may be determined by making the simple observation [24| that 

< 0^{t) >= TrF{0^{t)p^) = TrF{0^{t)TrA{p^+^)) = Tta+f{0^®I^P^^^) (14) 

where F and A correspond to the field and atomic variables, are the field operators and corresponds to the 
identity operator in the atomic space. We consider all expectation values to be normal ordered in the operators. This 
ensures that the spurious contribution of the vacuum is eliminated. The action of the creation and annihilation field 
operators on the atom+field states is given below : 

Ofcj |i; fc2; fci > ^0 for ^2 = 1 and V {i, fci}; 

= for fca = and V {i, fci}; (15) 

The expectation value of the annihilation operators ak',k~ 1,2 of the two modes is given by 

< ak{t) > = TrAMak{t)p'^+^) 

ki /c2 

X 04i^^(t)N;fc2;fci ><^;fc2;fcl|vI'><vI/| (16) 
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where O^^lifc'fc' ^ -^k', ®ik{ ® a^^ ®iki and ik2{ki) is the identity operator of mode ^2(^1) of the field. 

Identifying ^2 = 2, fci = 1 the expectation values simplify to 

<aiW> = E"HO«rW+^af(t)«r(t)=a}(tK*(t) 

i i 

<a,{t)> = Y.aUt)aT{t)+Y.^UtH'{t)=al{i)a^;{t)+al{t)a\'{t) (17) 

i i 

since the complex amplitudes {a\,a\,ci\,ci\^a\^a^} = due to RWA. The expectation values of the higher order 
correlation functions, specifically the normal ordered field correlation and the intensity correlation functions which 
are nonzero, are given by 

<a\{t)a,{t)> = Y.\a]{t)\^ = \a\{t)\^ 

i 

<al{t)a2{t)> = \al{t)\^ + \al{t)\^ 

i 

< alit)a\it)ai{t)a2it) > (18) 

We observe that < al{t)a\{t) > is the anomalous coherence [1^ or the cross correlation of the amplitudes of the 
two fields emitted by the atomic system. On the other hand, the expectation value of the Pauli operators for the two 
transitions is listed below: 

<a«(t)> = ^a](t)af (i)=a^(Oaf W 

7 

<a(^\t)> = 5]aI(0al'W=«^W"rW- (19) 
7 

The expectation values of the normal ordered second and the fourth order spin correlation operators are given by 

7 

< a^^\t)a^'Ht) > = 5:K(i)P = |a§(0P 

7 

7 

< a^_^\t)a^^\t)ai'\t)ai'\t) > = ^ \aUt)\' = \al{t)\' . (20) 

7 

From the above equation we see that < a'^^\t)a':^\t)a'L\t)a'L\t) > =< a':^\t)a^^\t) > and the equivalence between 
the field operators and the atomic transition operators imply the following set of identities: 

\al{t)f = 0; 4{t - r/c) oc a?(t); a^{t - r/c) oc a}(<). 

The first of the conditions is a manifestation of the Schmidt decomposition which forces to be of rank three. Note 
that the atomic operators are related to the amplitudes {a^, a^, 03} only while the emitted radiation is related to the 
amplitudes {a°, a}, } only. As already mentioned depending on the practical feasibility, one may choose observables 
from either the field space or the atomic space. 

C. Numerical results 

We now discuss the numerical results in three regimes of interest: a) the pure state regime fli < fl2, b) the Bell 
state regime, fii ~ fl2 and c) the mixed state regime f2i > f22. The reason for this nomenclature becomes clear in 
our analysis below. All the system parameters may be rendered dimensionless by scaling with the atomic lifetime 
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FIG. 2: (Color online) Time evolution of the populations for (a) pure state regime, Qi = 3.0, Q2 = 6.0 (b) the Bell state regime 
Qi = 6.0, fl2 ~ 6.0 and (c) the mixed state regime Qi = 6.0, ^2 = 3.0. The black dotted curve corresponds to Trp^ (purity). 
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FIG. 3: (Color online) Real part of the two photon density matrix for (a) the pure state regime, (c) the Bell state regime and 
(e) the mixed state regime and (b), (d), and (f) show the imaginary part of the density matrix in the respective regimes. The 
parameters are the same as in Fig. 2. 



7 ?» IMHz . In particular we choose the pure state regime corresponding to the driving field strengths fix = 3.0 and 
5^2 = 6.0. Our numerical simulations correspond to the two-photon resonant situation; therefore, Ai = A2 = 0. The 
upper most level |3) (which may correspond to 51)5/2) is assumed to be metastable, i.e., = 1.0 and r2 = 6.0. We 
further assume the system to be closed. 

Fig. 2 shows the populations in the |00), |01) and the |11) of the two photon state as functions of time. We have 
studied the time dependence of the populations in the three regimes mentioned above. Fig. 2a corresponds to the pure 
state regime wherein the population is negligible in |01) and is mostly in |00) with a relatively small admixture of 
Indeed it is clear from the graph that the purity, Tr{p'^), is quite close to one. The Bell state regime corresponds 
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to fii = 6.0, Q.2 = 6.0 and is illustrated in Fig. 2b. Here, the level [01) starts getting populated but remains small 
compared to the populations in |00) and the |11) states which become comparable to each other. The coherence 
between |00) and |11) is also quite large. In the mixed state regime (fii = 6.0, ft2 = 3.0) shown in Fig. 2c, we note 
that all the populations are comparable and the purity is much less than one. Note that the time evolution in all the 
regimes attains steady state in approximately an atomic lifetime after performing one Rabi cycle. This justifies the 
single photon approximation approach used here. 

The steady state two-photon density matrix corresponding to the three regimes is shown in Figs.3a-3f . The real 
part of the density matrix in the pure, the Bell and the mixed state regimes are respectively shown in Figs. 3a, 3c and 
3e. Likewise the imaginary part of the density matrix in the three regimes are shown in Fig. 3b, 3d and 3f. Note that 
the third column and the third row are zero everywhere, as implied by Schmidt decomposition and the dynamics. 
In the pure state regime (Figs. 3a, 3b), even though most of the population is in the ground state there is a nonzero 
population in the |11) state and a nonzero correlation in the Poo ii element. This is an almost pure state as is clear 
from Fig. 2a {Trp^ ss 1) and the pure state component is of the form ci|00) + C2III); ci 2> C2. The small mixedness 
is due to the nonzero population in the |01) state. In the Bell state regime (Figs. 3c, 3d), the correlation is strong and 
further the populations in the |00) and |11) are equal . The loss of purity is again due to the |01) state. In the mixed 
state regime (Figs.3e,3f), the steady populations in the three levels are comparable, and the correlations /Ooo;0i and 
Poi;ii compete with the correlation poo;ii- This regime provides the most general two photon mixed state . 

Having determined the state of the two-photon density matrix we now move on to study the entanglement properties 
of the two-mode emitted radiation. 



IV. CONTROL AND CHARACTERIZATION OF ENTANGLEMENT 



Pure state entanglement may be characterized by various measures such as the Von Neumann entropy of the 
reduced density matrix, concurrence and negativity. All these measures are equivalent in the sense that they are 
relative monotones with respect to each other. In contrast, it is well known that the entanglement of a mixed state 
(MSE) cannot be characterized by a single parameter [2^ . Generalizations of the above mentioned pure state measures 
in terms of a single parameter cease to be relative monotones of each other in the mixed state regime; this has been 
demonstrated explicitly in the case of concurrence and negativity [l^]. In a detailed and exhaustive analysis, Peters 
et. al. [2^, [2^ have analyzed the sensitivity of the state to a number of measures of entanglement for depolarized 
nonmaximally entangled states such as the Werner states and also the so called maximally entangled mixed states. 
They perform a relative comparison of six benchmarks of entanglement viz., fidelity, trace distance, tangle, linear 
entropy, concurrence and von Neumann entropy, and find an imbalance in the values of these quantities, even in 
situations where it is not expected, clearly signalling that none of these can by itself fully account for MSE. Finally, 
we mention quantum discord [30j which again displays the tension between such definitions of entanglement. 

We would like to point out here that the description of mixed state entanglement by Bhardwaj and Ravishankar 
(Tsj addresses these issues by characterizing MSE in terms of a probability density function (PDF) rather than a single 
number. This characterization does not invoke any new concept such as entanglement of formation or separability. 
Instead, it is based on the observation that a description in terms of mixed states is required when we have an ensemble 
of quantum systems each of which is in a pure state. A probabilistic description is, therefore, natural. This expectation 
is vindicated by the fact that while two states which are not equivalent under local operations {SU{2) ® SU{2)) can 
have the same value for concurrence, the probability densities are necessarily distinct. In fact, the state can be 
reconstructed almost uniquely (up to local transformations) by the density. We employ this characterization of MSE 
in this paper. We briefly describe their formalism: 



A. Probability density function for MSE 



The characterization of mixed state is accomplished in two steps. Firstly, the special case when the state is a 
projection of rank d is considered (For a two qubit system, dmax = 4), = ^^d, where the operator Prf projects a d 
dimensional subspace Hp^ of the two qubit Hilbert space ^4. We then observe that for all tp € Hd, — Ij which 

means that there is a uniform probability for the system to be in any state in the subspace "Hd- By the same token, 
for all states tp ^ "Hd, {'<P)pd — 0, implying that the probability for finding the system in the complementary subspace 
vanishes identically. Armed with this simple basis independent observation, define the probability density function 



JdVnA£^~£) 
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where dVn^ is the volume measure for the manifold Hp^. To repeat, the density so obtained is not an artifact of any 
expansion of p. We mention that the state may be reconstructed from the resultant PDF, which is unique — up to 
local transformations. 

In the case when p is not a projection, it is written as a weighted sum of projection operators P^, satisfying the 
nestedness condition 'Hp^_j C Hp^. Consider the expansion of p in its eigenbasis with its eigenvalues written in a 
non-increasing order, then, 

1=1 d=l ^ 

with the restriction that W4 = A4/A1. Employing this resolution, the PDF for any state can be written as 

4 

Pp(f)=^wrf7'p,(f). (23) 

d=l 

The PDF so defined is invariant under SU{2) SU{2) transformations. 

It has been shown by Bhardwaj and Ravishankar that it requires seven independent parameters which characterize 
the PDF in various subspaces as follows: 

Vp, ^ W4. (24) 

Recall that pd is the projection onto the d dimensional subspace. In what follows, for pure state entanglement, we 
shall employ the pure state concurrence, £ = 2\{'ip'^ where {if)'^) is the time reversed state of 



1. PDF for projective states 

Before we proceed to use the above characterization of mixed state entanglement, a brief description of the param- 
eters is in order. 

(i) Pure states: Consider first a state which is a one dimensional projection. Being pure, it has a sharp value of 
entanglement. The PDF is, accordingly, given by a Dirac delta: 

Vp, ^5{£-£^)- Pi = (25) 

(ii) Two dimensional projections p2 '■ The evaluation of PDF for a two dimensional space involves two parameters and 
is somewhat involved, given its rich structure. We briefly outline the idea. Consider a state l^*) G ^{{¥2) which may, 
in general, be written in terms of two orthonormal states |xi), IX2) as 

|vl/) ^ cos ^e'^/2|xi) + sin ^e'^^/^lxs) (26) 

where < 6 < n and Q < cj) < 2tt. It follows from elementary algebra that one can always construct a linear 
superposition which is separable, which can be conveniently be labelled |00), after a further local SU{2) (3 SU{2) 
transformation. Using this as a basis vector, its orthogonal state x|01) +y|10} -|-z|ll) may be so chosen, by employing 
the residual SO{2) ® 50(2) symmetry that preserves |00), such that x, y, z are real nonnegative. By the normalization 
condition, only two of them, say x, y are independent. In short any state which is a two dimensional projection is 
characterized by two invariant parameters under local transformation. 

The mixed state density matrix of such a state (two qubits A and B) may also be written in the form 

= {l + dA-VA + ^B-^B + ® Ob) ■ n) (27) 

Here it may be shown from the equivalence of the atomic and field operators that Pa and Fb are the polarizations 
due to the two transitions of the atomic system; in particular we identify A — >■ 1,5 — > 2. The invariant parameters 
X, y may be written in terms of the polarizations as 

x^ = Pl/{1 + ^l-Pl-Pl) 

= pya+^^-pi-pi) 

= (l-x2-y2) (28) 
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FIG. 4: (color online) Invariant parameters and correlation as functions of the ground state excitation energy Sli; ^2 ~ 3.0. 
The parameters wl, w2 and w3 correspond to the three weights w;, 'entl' = £i, 'ecusp' =£ cusp, 'emax'= iPmax are defined in 
eq.(24) and 'corr'= Pog ii 



where PX(b} = ^ ~ 4:det{p^'-^^) and = Trsip^^)- The PDF is shown to attain the form 

Here £ satisfies the bounds < f < U{e), C{9) < £ < where 14(6) = \zsm9 + (1 - cos6'(l - z'^)/2\ and C{9) = 
zsinfl — (1 — cos0(l — z^)/2|. Note that in the above equation the Haar measure for the volume is simply that of the 
two sphere S^: dVua ~ smOdOAcj). 

Three dimensional projections p^: The PDF for three dimensional projections have a more complicated volume 
measure, but a simpler structure, made possible by the fact that it is uniquely characterized by its dual, the orthogonal 
(pure) state \ip±). Let its entanglement be £±. The PDF for p^ is then given by 



Pp3 {£) = _^ cosh-i — (30) 



where £"> = max(£, £±). As an interesting aside, we note that concurrence and negativity vanish identically for any 
three dimensional projection. 

The four dimensional projection is simply the fully unpolarized system, which as a universal background. We 
do not reproduce the curve here since it has zero weight for our system since one of the eigenvalues of the photon 
density matrix, corresponding to the state |10), is zero. There is a corresponding reduction in the number of invariant 
parameters that characterize entanglement. 

An equivalent description to PDF is given by the cumulative distribution function (CDF) defined by 

J-p(f) = f d£'Vp{£'), (31) 
Jo 

We would like to mention here that we have evaluated both the probabiUty density function (PDF) and the cumulative 
distribution function (CDF) for the entanglement as we find that the latter displays the dynamical features more 
clearly. 



B. Numerical Results and Discussion 



We have evaluated the steady state values of the six invariant parameters characterizing the mixed state entan- 
glement distribution, viz., the three weights wi,W2,W3, the entanglement of the pure state component £i, and the 
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FIG. 5: (Color online) Probability distribution and the cumulative distribution of entanglement corresponding to one, two 
and three dimensional subspaces. for the a) pure state, b) the Bell state and c) the mixed state regimes. Legend:( d) 
probability distribution function for one dimensional sub-space and (e) the corresponding cumulative distribution, (f) PDF for 
two-dimensional sub-space and (g) its CDF, (h) PDF corresponding to the three dimensional subspace and (i) its corresponding 
CDF. 
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two entanglement parameters £cusp,£max corresponding to the two dimensional projection. The parameters W4, £±^ 
vanish identically. The parameters £maxj £cusp a-rc determined from p2 (see Eq. 28) by 

£7nax = Xy + \J z"^ +x2y2 

£cusp ^ /£max- ('^^) 

The three interesting regimes, VLi < Q.2, ^1 ~ ^^2, and Q,i > Q2, create the two photon mode system in different 
states. In the first regime the state has a high degree of purity, but little entanglement. In the second regime, most 
appropriate for fine tuning, the system loses purity to an extent, but gains in entanglement, as indicated in Figs. 2 
and 3. The last regime destroys both purity and entanglement. Fig. 4 shows the variation of these parameters as 
functions of the Rabi frequency of the ground state excitation ili, with 0,2 = 37, where all the three regimes show up. 
But beyond this mildly interesting qualitative feature, it is clear that in much of the range, the individual parameters 
by themselves give little information on the nature of entanglement reflecting the fact that no single parameter 
can completely capture all the aspects of mixed state entanglement. In contrast, as we show below, the probabiity 
distribution characterized by these invariant parameters sheds light on not only the nature of mixedness of the state, 
but also the nature of entanglement without any ambiguity. This description seems to overcome the difficulties raised 
in references [2^ [2^. In Fig. 5 we show the PDF and the CDF of entanglement corresponding to each of one, two 
and three dimensional subspaces for the three regimes of interest. It is clear from Fig. 5a which corresponds to the 
pure state regime, the probability distribution has a dominant contribution of « 95% due to the one dimensional 
subspace. The corresponding state is not fully entangled, with £ « 0.7 . Recall that for a Bell state £ = 1. As to 
the precise nature of the state, one can infer from Fig. 3a that the state must be a superposition of the |00) and |11) 
states with an excess of population in the |00) state. This is further substantiated by the strong correlation as shown 
by the poo;ii element. The deviation from purity is due to the unequal populations in the |00) and the |11) state. 
The two and three dimensional contributions arc negligibly small in this case. In the Bell state regime. Fig. 5b, there 
is a dominant contribution from the pure state even though the probability decreases to « 80%. Note however, that 
the state in this regime is closer to a Bell state since £ is close to one. There is a small contribution from the three 
dimensional subspace. In the mixed state regime Fig. 5c there is a comparable contribution from all the subspaces. 
The pure state probability has however reduced to less than 50% and also £ « 0.65. The nature of mixedness in this 
regime may be understood by looking at Fig. 3c and Fig. 5c. Even though the populations in |00) and |11) are nearly 
equal, the correlation /9oo;ii is much weaker. On the other hand the correlation poo;Oi is stronger implying that the 
superposition of the |00) and |01) is competing. Since both the correlations arc comparable the states belong to either 
two or three dimensional subspace. 

The equivalence of the photon and atom density matrix suggests that the correlation /3oo ii(0 is proportional to 
the atomic coherence pf^{t — r/c). Note that maximum entanglement of the two photon occurs when this correlation 
and hence the atomic coherence is strong. This clearly demonstrates the fact that the quantum correlations of the 
atomic system manifest as strong correlations in the emitted radiation as well. In the mixed state regime the single 
and two photon coherences compete and results in a weaker entanglement of the two photons. Thus the control of 
entanglement/purity is achieved by controlling the atomic coherences. 

To conclude, we have shown that it is possible to determine the two photon density matrix, in its Fock space, by 
a combination of observations that include measurements on the atomic system that generate the photons and the 
radiation itself. By a proper choice of the two driving field strengths, the radiation emitted can either be prepared as 
a pure state, but with little entanglement, or with a reasonable overlap with the Bell state, but with some mixedness. 
Both the regimes are of interest. The former is eminently suited for the preparation of fiducial states as in applications 
such as the Deutsch algorithm. The latter regime, where highly entangled states are produced is useful when initial 
entanglement is required, as in the case of teleportation. The third regime, where fti ^ ^2 seems unsuited for 
applications; on the other hand, it provides mixed states exhibiting a rich variety of entanglement which cannot be 
captured in terms of a simple parameter. It should be of interest to extend this study to multi photon production 
and its characterization in terms of purity and entanglement. 
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